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Picea koyamae is a species which is distributed scarcely and is threatened to extinction, though the details 
are not well known. We investigated the geographical distribution and genetic variation of this species. 
Seven habitats over the Yatsugatake mountains and Akaishi mountains were confirmed in this study, and 
it was estimated that in these areas combined there were less than 1,000 living mature trees. Genetic 
diversity within and among 4 populations was analyzed using allozymes. The mean genetic variation 
within populations was in the range of general long-lived tree species. However, genetic variation in the 
Karamatsu-sawa population at Yatsugatake mountains showed a rather low value. It is conceivable that 
this population has experienced a bottleneck effect. The genetic differentiation among 4 populations was 
found to be at an intermediate level between continuously distributing conifers and discontinuously dis¬ 
tributing ones. In this study, “P. shirasawae' ’ was suggested to be merely a morphological variation with¬ 
in a population of P. koyamae, and was treated as P. koyamae in our results. 

Key words: Akaishi mountains, allozyme, genetic diversity, habitat, Picea koyamae, Picea shira¬ 
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Picea koyamae Shiras. (Yatsugatake-tohi in 
Japanese) is a tall conifer tree species which is dis¬ 
tributed in Honshu Island, Japan. In 1911 it was dis¬ 
covered in the Yatsugatake mountains area 
(Shirasawa & Koyama 1913). Since then, P. koya¬ 
mae has been believed to be distributed only on 
Yatsugatake mountains (Farjon 1990, Yamazaki 
1995). Currently, there are less than 500 living 
mature trees in the Yatsugatake mountains area, 
and is listed as “threatened with extinction” in a 


report on Red Data Plants (Farjon & Page 1999, 
Environmental Agency of Japan 2000). 

The genus Picea is one of the most important 
tree groups in Japanese flora to come out of the 
Last Glacial Age, and fossils classified into P koya¬ 
mae or “P. shirasawae Hayashi” have been discov¬ 
ered in very wide area on Honshu Island (Minaki 
1987). In addition, more than 10 genus Picea 
species, including now-extinct ones, have been 
reported based on fossil records in Honshu (Suzuki 
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1991, LePage 2001), though the details of their 
phylogenetic relationships are unknown. 

There are several problems with the classifi¬ 
cation of Picea koyamae and its sibling species. 
Around the Japan Sea, six or seven species classified 
into Sect. Picea have been reported by Farjon 
(1990), Shimizu (1992) and Yamazaki (1995) (Table 
1, Fig. 1). Some botanists have pointed out that P 
koyamae, P. koraiensis, and P. obovata are very 
similar (Wilson 1916, Wright 1955). In particular, 
there was a case that P koraiensis was included 
with P. koyamae as one species (Fu et al. 1978), 
though no detailed examination had been made. 


El 20" El 50" 



Fig. 1. Distribution of Sect. Picea species around the Japan Sea. 
The species names are followed by description of Farjon 
(1990). 


On the other hand, “ Picea shirasawae ” (Hime- 
matsuhada in Japanese) has been sometimes treated 
as an independent species (Hayashi 1960b, Yama¬ 
zaki 1995). However, many botanists have classified 
it as a variety of P. koyamae (Shimizu 1989) or P. 
alcoquiana (Faijon 1990). “Picea. shirasawae ” was 
discovered in the same geographical area in which 
P. koyamae was discovered, and the discoveries of 
both were published in the same report (Shirasawa 
& Koyama 1913). “Picea. shirasawae ” resembles P. 
koyamae, and has been reported to be morphologi¬ 


cally distinguished from P koyamae by longer cone 
length (8-12 cm), longer leaf length (13 - 18 mm) 
and branchlets color (orange-brown in pubescence) 
(Hayashi 1960a). However, Shimizu (1992) pointed 
out that there was not a difference other than the 
form and size of cone between those two taxa. In 
addition, Kobayashi et al. (2000) reported a total of 
about 2500 bases of nucleotide sequences in four 
regions of chloroplast DNA for one P. koyamae 
individual and one “P. shirasawae ” individual, and 
their sequences were identical with each other 
though they showed some difference from the 
sequences of other four Japanese Sect. Picea 
species. This might support the classification of 
Shimizu (1992). 

In this study, we confirmed the present geo¬ 
graphical distribution of Picea koyamae and ana¬ 
lyzed their genetic variation using allozymes 
(enzyme polymorphisms) without distinguishing 
“P. shirasawae ” as a separate species according to 
the classfication by Shimizu (1989, 1992). 

Materials and Methods 

Plant materials 

Based on descriptions in past reports (Hayashi 
1960a, Yamazaki 1965, Katsuki & Seido 1999), 
we considered Picea koyamae to be geographically 
distributed in the Yatsugatake mountains and 
Akaishi mountains areas. Then, we confirmed their 
current habitats by means of field survey and direct 
examination to identify the species. 

After this survey, we studied the genetic vari¬ 
ation in allozymes of the four populations (Fuki- 
sawa and Karamatsu-sawa populations in Yatsuga¬ 
take mountains; Ohdaira and Tenshi-iwa popula¬ 
tions in Akaishi mountains) from which we had 
obtained a sufficient number of samples of Picea 
koyamae individuals (Fig. 2). 

For allozyme analysis, we collected young 
branches of about 30 individuals from each of the 
four populations in June, 2000. Collected branches 
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were stored at -80 °C until the enzymes were extract¬ 
ed. 

Allozyme analysis 

Procedures for enzyme extraction and polyacry¬ 
lamide gel electrophoresis followed those described 
in Tsumura et al. (1990) and Tsumura & Ohba 
(1993), except that 100 mg of young leaf tissue 
from each tree and 80 mg of polyvinyl-polypyrroli- 
done were ground to a paste with a mortar and pes¬ 
tle in 1 ml of extraction buffer (100 mM Tris-HCl 
pH 7.5, 23.4% glycerol, 0.6% tween 80, 1.2 mM 
dithiothreitol, 0.6 mM NAD, 0.5 mM NADP, 0.5% 
2-mercaptoethanol, and 0.08% albumin). 

Eight enzyme systems were employed: 
diaphorase ( Dia , EC 1.6.4.3.), fumarate hydratase 
(Fm, EC 4.2.1.2.), glycerol-2-dehydrogenase ( G2dh , 
EC 1.1.1.29.), glutamate dehydrogenase (Gdh, EC 
1.4.1.2.), glutamate oxaloacetate transaminase (Got, 
EC 2.6.1.1.), phosphoglucomutase ( Pgm , EC 
2.7.5.1.), shikimate dehydrogenase (Shdh, EC 
1.1.1.25.), and UDP glucose pyrophosphorylase 
( Ugpp, EC 2.7.7.9.). The staining method for Ugpp 
is described in Harris & Hopkinson (1976), and 
those for the other enzyme systems follow the meth¬ 
ods set out by Tsumura et al. (1990) and Tsumura & 
Ohba (1993). 

Statistical analysis 
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Levels of genetic variation were calculated for each 
of the populations according to a number of para¬ 
meters including: the proportion of polymorphic 
loci (P) at the 95% criterion, the average number of 
alleles per locus (Na), the effective number of alle¬ 
les per locus (Ne; Kimura & Crow 1964), the 
observed heterozygosity (Ho), and the expected 
unbiased heterozygosity under Hardy-Weinberg 
equilibrium (He; Nei 1978, Nei & Roychoudhury 
1974). 

F-statistics (Fis) (Wright 1951) that represent 
the deviation of genotypic proportions from the 
Hardy-Weinberg expectation were calculated for 
polymorphic loci in each population. Statistical 
analysis of deviations from such expectations was 
made using the % 2 test (Li & Horvitz 1953). 

Gene differentiation within and among popu¬ 
lations were calculated according to the following 
genetic diversity statistics (Nei 1973, Nei & Chesser 
1983): the total population gene diversity, Ht; the 
average gene diversities within populations, Hs; 
and the relative extent of gene differentiation among 
populations, Gst. 

Nei’s unbiased standard genetic distances were 
calculated for all population pairs (Nei 1972, Nei & 
Roychoudhury 1974), and clusters were analyzed 
using the unweighted pair-group method with arith¬ 
metic averages (UPGMA; Sneath & Sokal 1973). 

The above calculations were carried out using 
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Table 2. The confirmed populations of Picea. koyamae in this study. 


Mountains 

Habitat 

Location 

Altitude 

(m) 

Population 
size * 

Yatsugatake 

Fuki-sawa 

Nishidake-305 §§ , Fujimi, Nagano 

1800 

30 

Karamatsu-sawa 

Nishidake-310 §§ , Fujimi, Nagano 

1650 

100 


Ohdaira 

Nirasaki-66 §§§ , Hakushu, Yamanashi 

1600 

10 


Kitazawa-toge 

Miwa-278 §§ , Hase, Nagano 

1500 

<10 

Akaishi 

Mikobuchi 

Inasato-78 §s , Hase, Nagano 

1450 

<10 


Koma-sawa 

Inasato-85* 5 , Hase, Nagano 

1450 

<10 


Tenshi-iwa 

Ohkawara-76" 5 , Ohshika, Nagano 

1750 

50 


* Roughly estimated number of mature trees (bigger than about 20cm DBH). 
Forest compartment assigned by the National Forest Agency. 

Forest compartment assigned by Yamanashi Prefecture. 
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GENESTRUT (Constantine et al. 1994) and DIS¬ 
PAN (Ota 1993) computer software. 

Results 

Present geographical distribution 
In our survey, a total of about 220 mature trees of 
Picea koyamae at 7 habitats were confirmed (Fig. 2, 
Table 2). 

In the Yatsugatake mountains area, Fuki-sawa 
was thought to be the type locality of the Picea 
koyamae collected by Mr. Koyama in 1911 (Wilson 
1916), but this forest suffered heavy damage by 
the Isewan typhoon in 1959. At present, over 1,000 
young trees and about 30 mature trees (bigger than 
about 20cm DBH) with cones can be observed in 
Fuki-sawa. Karamatsu-sawa, where there are about 
100 mature trees of P koyamae, is the most well 
known habitat (Hayashi 1960a). The National Forest 
Agency registered it as a gene resource preservation 
forest in 1987. In this population, the diameter dis¬ 


tribution of mature trees shows a peak at 30-40cm. 
The age of a big tree (DBH = 50.0cm; H = 22.3m) 
damaged in 2001 by typhoon was estimated to be 
103 years. 

In the Akaishi mountains area, Ohdaira is a 
locality that has been confirmed recently as a habi¬ 
tat of Picea koyamae and P. maximowiczii (Katsuki 
& Seido 1999), and about 30 individuals of P koya¬ 
mae with a height over 1.3m were observed in a 
small area (14ha). The number of mature trees 
seems to be about 10. Kitazawa-toge, Mikobuchi, 
and Koma-sawa are the localities reported by 
Hayashi (1960a). However, the numbers of samples 
of P. koyamae in those habitats were insufficient for 
the present allozyme analysis, since the number of 
mature trees in each of those localities appeared 
to be less than 10. Tenshi-iwa was reported by 
Yamazaki (1965) to be the southern limit locality of 
“E shirasawae , \ Though the details of this popula¬ 
tion were not reported in the previous paper, our 
field survey clarified that there are more than 50 
mature trees of P. koyamae in this locality. 


El80° 00' 


El80° 30' 


N36° 00 1 


N35° 30’ 


Nagano 


Pref. 


Yatsugatake mountains 


Fuki-sawax 
Karamatsu-sawa A 


X Ohdaira 


Kitazawa-toge q 


Akaishi mountains 



Fig. 2. Reported distribution and confirmed habitats of Picea 
koyamae. Hatched area, geographical distribution described 
in Hayashi 1960a, Katsuki & Seido 1999 and Yamazaki 
1965. X mark, habitat in which samples were collected. 
Circle, observed habitat. 


Genetic variation in 4 populations at two moun¬ 
tainous areas 

The sampled individuals for allozyme analysis were 
mainly mature trees in Fuki-sawa, Karamatsu-sawa 
and Tenshi-iwa and the mixture of mature trees 
and young trees in Ohdaira. Ten putative loci were 
detected on eight enzyme systems: two loci each for 
Dia and Got and one locus for the others. At the 
95% criterion, four loci {Dia-2, G2d, Pgm, and 
Shdh) were polymorphic; one locus ( Dia-1 ) was 
classified as monomorphic with very little varia¬ 
tion; and the other five loci (Fm, Gdh, Got-1, Got-2, 
and Ugpp ) showed no genetic variation at all (Table 
3). Table 4 shows a summary of the degree of genet¬ 
ic variation. It can be seen that three populations 
(Fuki-sawa, Ohdaira and Tenshi-iwa) showed 
approximately the same degree of genetic varia¬ 
tion (He=0.185-0.196). The Karamatsu-sawa pop¬ 
ulation appeared to hold slightly lower variation 
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(He=0.136) than those three populations, though 
the difference was not statistically significant. 

The estimated values of Fis for polymorphic 
loci are summarized in Table 5. A significant Fis 
value is seen only for the Pgm locus in the Ohdaira 
population; other Fis values were not significant. 

Table 6 shows the levels of genetic differenti¬ 
ation among the 4 populations. The average Gst 
value over all loci was 0.067. 

Genetic distances among the 4 populations 
were between 0.001 (Ohdaira vs. Tenshi-iwa) and 
0.045 (Fuki-sawa vs. Karamatsu-sawa) (Table 7). 
Based on standard genetic distances, the dendrogram 
by the UPGMA method is shown in Fig. 3. The 
bootstrap probability suggested combining the 
Ohdaira and Tenshi-iwa populations as the most 
closely related at 69%. 

Discussion 

In this study, 7 habitats of Picea koyamae were 
confirmed (Table 2), but the number of trees in 
Akaishi mountains area were shown to be underes¬ 
timated, since our survey was not conducted in all 
areas of Akaishi mountains. Furthermore, in some 
cases, a stand of P. koyamae individuals occurred on 
a steep slope, and we could not reach such a habitat. 
Note et al. (1999) reported many P koyamae indi¬ 
viduals standing on steep slopes or near cliffs on the 
other side of the Mibu River at Mikobuchi by view¬ 
ing them with binoculars, but it was not possible to 
approach these trees and identify their species direct¬ 
ly. Table 2 includes only the number of mature 
trees we could identify directly by touching them. 
Consequently, it is supposed that in the area of the 
Akaishi mountains there are more than 500 mature 
trees that can produce cones. This estimation sup¬ 
ported inclusion of this species as one “vulnerable to 
extinction” on the red list by the Environmental 
Agency of Japan (2000). However, population size 
and density are not sufficient for regeneration in 
habitats where mature trees stand alone, and it is 


feared that small populations will disappear in the 
future. 

With regard to genetic variation within the 3 
populations (Fuki-sawa, Ohdaira and Tenshi-iwa) of 
Picea koyamae, shown in Table 4, estimated values 
of Ne (1.23-1.25) and He (0.185-0.196) in this 
study were not clearly different from the summa¬ 
rized values (Ne=1.32, He=0.173) of 195 long- 
lived tree species (Hamrick & Godt 1992). In com¬ 
parison with Japanese conifer species, the genetic 
variation of P koyamae is nearly the same as those 
of Cryptomeria japonica (Ne=1.33, He=0.189; 
Tomaru et al. 1994), Chamaecyparis obtusa 
(Ne=1.33, He=0.193; Uchida et al. 1997), and Pinus 
pumila (Ne=1.38, He=0.225; Tani et al. 1996), and 
greater than that of Abies mariesii (Ne=1.08, 
He=0.054; Suyama et al. 1997). On the other hand, 
variation of the Karamatsu-sawa population 
(Ne=1.16, He=0.136) was rather lower than those of 
the other populations investigated in this study. In 
this population, diameter distribution of mature 
trees shows one peak, and the forest age was esti¬ 
mated to be about 100 years old. An investigation 
conducted in 1911 showed no old growth at that 
time (Shirasawa & Koyama 1913). From those evi¬ 
dences, this population might have regenerated at 
once, about 100 years ago. It is conceivable that the 
low variation within this population would have 
been caused by a bottleneck effect through a limit¬ 
ed number of mother trees. 

With regard to genetic differentiation among 
the 4 populations of Picea koyamae shown in Table 
6, the estimated Gst value (0.067) in this study 
was slightly lower than the mean estimated value 
(0.084) summarized for 195 long-lived tree species 
(Hamrick & Godt 1992). Gst values of Japanese 
conifers were estimated to be 0.034 for Cryptomeria 
japonica (Tomaru et al. 1994), 0.030 for Chamae¬ 
cyparis obtusa (Uchida et al 1997), 0.170 for Pinus 
pumila (Tani et al. 1996), and 0.144 for Abies 
mariesii (Suyama et al. 1997). The former two 
species have wide (over all of Japan except 
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Table 3. Allele frequencies of 10 allozymc loci in 4 populations. 


Yatsugateke Akaishi 


Locus - 

Allele 

Fuki-sawa 

( 3 1§§§) 

Karamatsu-sawa 

(25) 

Ohdaira 

(32) 

Tenshi-iwa 

(31) 

Dia-P 


0.97 

1.00 

0.97 

0.97 


b 

0.03 

0.00 

0.03 

0.03 

Dia-2 

a 

0.60 

0.08 

0.30 

0.47 


b 

0.40 

0.92 

0.70 

0.53 

Fm 

a 

1.00 

1.00 

1.00 

1.00 

G2d 

a 

0.42 

0.38 

0.36 

0.45 


b 

0.58 

0.62 

0.64 

0.55 

Got-1 

a 

1.00 

1.00 

1.00 

1.00 

Got-2 

a 

1.00 

1.00 

1.00 

1.00 

Pgm 

a 

0.48 

0.70 

0.69 

0.68 


b 

0.52 

0.30 

0.31 

0.32 

Shdh 

a 

0.27 

0.18 

0.34 

0.32 


b 

0.73 

0.82 

0.66 

0.68 

Ugpp 

a 

1.00 

1.00 

1.00 

1.00 


f The locus with -1 in its designation represents a form moving electrophoretically 
faster than that with -2. 

■' 5 Allele a stands for the allele moving faster than allele b in electrophoresis. 

Sample number 


Table 4. Genetic variation within 4 populations. 


Mountains 

Population 

P(%) 

Na 

Ne 

Ho 

He(SE) 

Yatsugatake 

Fuki-sawa 

40 

1.50 

1.25 

0.213 

0.196(0.076) 


Karamatsu-sawa 

40 

1.40 

1.16 

0.132 

0.136(0.062) 

Akaishi 

Ohdaira 

40 

1.50 

1.23 

0.150 

0.185(0.072) 


Tenshi-iwa 

40 

1.50 

1.24 

0.184 

0.196(0.076) 


P, percentage of polymorphic loci at a 95% criterion; Na, number of alleles per locus; Ne, effec¬ 
tive number of alleles per locus; Ho, average observed heterozygosity; He, average expected 
heterozygosity. 


Table 5. Fixation indices (Fis) within 4 populations for polymorphic loci. 


Mountains 

Population 

Dia-2 

G2d 

Pgm 

Shdh 

Yatsugatake 

Fuki- sawa 

-0.121 

-0.173 

0.110 

-0.196 


Karamatsu-sawa 

-0.065 

0.085 

0.160 

-0.195 

Akaishi 

Ohdaira 

0.042 

-0.136 

0.570** 

0.318 


Tenshi-iwa 

0.044 

0.103 

0.129 

-0.017 


** PO.Ol. 
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Table 6. Levels of genetic differentiation among 4 populations 
based on 4 polymorphic loci and on all loci including 
those that were polymorphic and monomorphic. 


Locus 

Ht 

Hs 

Gst 

Dia-2 

0.450 

0.353 

0.216 

G2d 

0.481 

0.478 

0.005 

Pgm 

0.462 

0.447 

0.034 

Shdh 

0.403 

0.395 

0.020 

Mean 5 

0.449 

0.418 

0.069 

All 

0.184 

0.172 

0.067 


4 Mean on 4 polymorphic loci. 


Hokkaido) and continuous distribution areas, the 
latter two species’ distribution areas are restricted to 
central and northern Japan and are discontinuous. 
The value of 0.067 in P. koyamae is larger than 
those of the first two species and smaller than those 
of the latter two species. The Gst value of P. koya¬ 
mae might have been affected by its wide-area con- 


Vol. 55 

tinuous distribution on Honshu Island in the Last 
Glacial Age. 

The dendrogram (Fig. 3) based on standard 
genetic distances suggests that Ohdaira and Tenshi- 
iwa populations in Akaishi mountains are very 
closely related, and that the Karamatsu-sawa popu¬ 
lation is remote from the other three populations. 
The two populations in Akaishi mountains might be 
continuous through some unknown populations in 
areas which our field survey did not cover. The 
genetic remoteness of Karamatsu-sawa from the 
other populations could be explained by the proba¬ 
ble bottleneck effect. 

The Karamatsu-sawa population is preserved as 
a forest reserve, but its genetic variation is low and 
cannot represent that of the species Picea koya¬ 
mae. Therefore, for the purpose of preserving the 
genetic variation across the species, not only 
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Table 7. Standard genetic distances for all pairs from each of 4 populations. 

Fuki-sawa Karamatsu-sawa Ohdaira 

Karamatsu-sawa 0.045 (0.042) 

Ohdaira 0.014(0.012) 0.010(0.009) 

Tenshi-iwa 0.003 (0.004) 0.024 (0.025) 0.001 (0.003) 

The figures in parentheses stand for standard errors of genetic distances. 



I 


i 


0.010 0.005 

Standard genetic distance 


0 


Fig. 3. Dendrogram from UPGMA analysis based on pairwise standard genetic distances. The figure in the phylogenetic tree presents 
bootstrap probability based on 1,000 replicates. 
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Karamatsu-sawa but also other populations such 
as Fuki-sawa or Tenshi-iwa are important. It is 
desirable that care will be taken to preserve the 
populations in these other habitats, too. Fortunately, 
the Ohdaira population is scheduled to be regis¬ 
tered as a preserved forest by Yamanashi Prefecture. 

“Picea shirasawae ” was included as P. koya¬ 
mae in this study. There was no clear genetic evi¬ 
dence to suggest the existence of two independent 
taxa in our results, as follows. If P koyamae and “ P 
shirasawae ” are independent taxa and had not mated 
with each other in any of our samples, the Fis cal¬ 
culated on the assumption of a single random mat¬ 
ing population should have a positive value accord¬ 
ing to “Wahlund’s Principle” (Crow & Kimura 
1970). Flowever, table 5 shows that there was no 
clear result of significantly positive Fis, especially in 
the Fuki-sawa and Karamatsu-sawa populations 
which were reported as mixed stands of P koyamae 
and “P. shirasawae ” by Shirasawa & Koyama 
(1913) and Hayashi (1960a). According to the 
National Forest Agency, 40 individuals of P. koya¬ 
mae and 31 individuals of “P. shirasawae ” were 
registered from Karamatsu-sawa population in 1965. 
This evidence would suggest that “P. shirasawae ” is 
not an independent taxon, but rather, represents a 
morphological variation within P koyamae. 

We cordially thank T. Maeda, Y. Tsumura, K. Yoshimura, 
Y. Nishiyama, and K. Shimada (Forestry and Forest 
Products Research Institute) for help at sampling and 
examination; S. Shiraisi (Kyushu Univ.), T. Tanimoto 
(Utsunomiya Univ.), M. Suzuki (Tohoku Univ.), and M. 
Matsumoto (Chiba Univ.) for useful comments; students 
of Forest Science Department of Tokyo Univ. of 
Agriculture for help at sampling; and the Nanshin District 
Forest Office, Chubu Regional Forest Office, and National 
Forestry Agency for sampling permission. 
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